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This paper considers Gaussian flows multiplexed in a queueing 
network. A single node being a useful but often incomplete setting, 
we examine more advanced models. We focus on a (two-node) tan- 
dem queue, fed by a large number of Gaussian inputs. With service 
rates and buffer sizes at both nodes scaled appropriately, Schilder's 
sample-path large-deviations theorem can be applied to calculate the 
asymptotics of the overflow probability of the second queue. More 
specifically, we derive a lower bound on the exponential decay rate 
of this overflow probability and present an explicit condition for the 
lower bound to match the exact decay rate. Examples show that this 
condition holds for a broad range of frequently used Gaussian inputs. 
The last part of the paper concentrates on a model for a single node, 
equipped with a priority scheduling policy. We show that the analy- 
sis of the tandem queue directly carries over to this priority queueing 
system. 

1. Introduction. Traffic engineering in communication networks greatly 
benefits from models that are capable of accurately describing and predict- 
ing the performance of the system. This modeling is a challenging task, as 
a broad variety of traffic types are multiplexed in the network, with each 
of them having its specific (stochastic) characteristics. A commonly used 
modeling step is to represent the network nodes as queues, and to use queue- 
ing theory to analyze the performance (in terms of loss, delay, throughput, 
etc.) of the nodes. For the single queue operating under the first-in-first-out 
(FIFO) discipline, even for advanced traffic models detailed analyses are 
available. Evidently this single-node FIFO model gives valuable insights, 
but is an oversimplification of reality We mention two serious limitations. 
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First, traffic streams usually traverse concatenations of hops (rather than 
just a single node). Second, it is envisaged that the service at these hops 
distinguishes between several traffic classes (by using priority mechanisms, 
or the more advanced generalized processor sharing discipline) ; compare the 
Differentiated Services (diffserv) approach proposed by the Internet Engi- 
neering Task Force [17]. This motivates the recent interest in performance 
evaluation for these more complex queueing models. 

As indicated above, each type of traffic has its own stochastic properties, 
often summarized by the correlation structure. Traditional traffic models al- 
low only short-range dependent traffic processes, such as Markov-modulated 
Poisson processes or exponential on-off sources, in which correlations de- 
cay relatively quickly. Traffic measurements in the 1990s, however, showed 
that in various situations long-range dependent traffic models are more ap- 
propriate. This explains the popularity of Gaussian models, as they cover 
both short-range (cf. Ornstein-Uhlenbeck) and long-range dependent mod- 
els (e.g., fractional Brownian motion, see [19]). Another complicating issue 
is the fact that network traffic is usually influenced by feedback loops (think 
of TCP), which control how the user's traffic supply is transmitted into the 
network. Kilpi and Norros [18], however, argue that (nonfeedback) Gaussian 
traffic models are justified as long as the aggregation is sufficiently large 
(both in time and number of flows), due to central limit type of arguments. 

This paper concentrates on the evaluation of tail asymptotics in queueing 
systems that are more advanced than a single FIFO node. More specifically, 
we examine in detail tandem queues (particularly the second queue) and 
priority queues (particularly the low-priority queue); it turns out that the 
analysis of the tandem queue essentially carries over to the priority system. 
Our paper is meant as a first step towards the analysis of networks with gen- 
eral topology, with nodes operating under advanced scheduling disciplines 
such as Generalized Processor Sharing (GPS). 

In the tandem model we assume that n i.i.d. Gaussian sources feed into 
the queueing system, where the (deterministic) service rates of the queues 
as well as the buffer thresholds are scaled by n, too. We now let n go to 
infinity; the resulting framework is often referred to as the many-sources 
scaling, as was introduced in [32]. 

A vast body of results exists for single FIFO queues under the many- 
sources scaling. Most notably, under very mild conditions on the source 
behavior, it is possible to calculate the exponential decay of the probability 
p n (b, c) that the queue (fed by n sources, and emptied at a deterministic rate 
nc) exceeds level nb. Early references in this large-deviations framework are 
the logarithmic asymptotics found in, for example, [7] and [8]. We remark 
that exact asymptotics for Gaussian inputs were recently found by Debicki 
and Mandjes [9]. For Gaussian sources the logarithmic asymptotics of [7] 
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read 



(1) 



lim — logp n (b, c) = — inf 

n— >oo n t>0 



(b+(c-/x)t) 2 
2v(t) 



where \x is the mean input rate per source, and v(t) is the variance of the 
amount of traffic generated by a single source in a time interval of length t. 
The goal of the present paper is to find expressions similar to (1) for tandem 
and priority queues. 

Our work fits in the framework of a series of articles by Mannersalo and 
Norros [1, 23, 24, 25]. These papers examine queues with Gaussian sources, 
such as the single-node FIFO queue, but also priority queues and queues 
with GPS scheduling. For the latter types of queues, they derive heuristics 
for the decay rate of the overflow probabilities. The present paper shows 
that, for priority queues, these heuristics are typically close, but that there 
is a gap with the exact outcome. For both the tandem and priority queue 
a lower bound on the decay rate of the overflow probability is derived. In 
addition, we present an explicit condition under which this lower bound 
matches the exact value of the decay rate. Notice that lower bounds of the 
decay rate are usually of practical interest, as typically the network has to 
be designed such that overflow is sufficiently rare. 

Our analysis exploits the above-mentioned similarity between priority and 
tandem queues. The techniques applied stem from large-deviations theory, 
particularly sample-path large deviations, based on (the generalized version 
of) Schilder's theorem. We mention that for priority systems in discrete 
time, different bounds were found by Wischik [33]; we will comment on the 
relation with our results later. 

The paper is organized as follows. Section 2 introduces the tandem model, 
and presents preliminaries on (sample-path) large deviations. Section 3 an- 
alyzes the decay rate of the overflow probability of the second queue in a 
tandem system. This analysis is illustrated in Section 4 by a number of 
(analytical and numerical) examples. Section 5 studies the priority system, 
addressing the decay rate of the overflow probability in the low-priority 
queue. 

2. Model and preliminaries. This section introduces the tandem model 
that is analyzed in Section 3. In addition, we present preliminaries on large- 
deviations theory and the many-sources scaling. 

2.1. Tandem model. Consider a two-queue tandem model, with (deter- 
ministic) service rate nc\ for the first queue and nc2 for the second queue. 
We assume that c\ > C2 , in order to exclude the trivial case where the buffer 
of the second queue cannot build up. 
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We consider n sources (whose characteristics are specified in Section 2.2) 
that feed into the first queue. Traffic of these sources that has been served 
at the first queue immediately flows into the second queue — we assume no 
additional sources to feed the second queue. We are interested in the steady- 
state probability of the buffer content of the second queue Q2, n exceeding a 
certain threshold nb, b > 0, when the number of sources gets large, or, more 
specifically, its logarithmic asymptotics: 



Note that we assume the buffer sizes of both queues to be infinite. We remark 
that it is not a priori clear that the limit in (2) exists; its existence is a result 
of our study (Theorem 3.1). 

2.2. Gaussian sources. Let Ai(-) denote i.i.d. centered Gaussian pro- 
cesses with continuous sample paths and stationary increments, and Ai(0) 
0. for i = 1, .. . ,n. Then, for s <t, we interpret Ai(s,t) := Ai(t) — Ai(s) as 
the amount of traffic generated by the ith source in (s,i\. Denote by A(s,t) 
the generic random variable corresponding to a single source. The Gaussian 
sources are characterized by their variance function v(-) (which is necessarily 
continuous); for s <t, Var A(s,t) = v(t — s). 

Although in this setup the Gaussian processes are centered, our analysis is 
capable of handling the situation in which the sources have a positive mean 
traffic rate fi (smaller than both c\ and C2, to guarantee stability). This 
is due to the fact that the results for centered sources can be translated 
immediately into results for noncentered sources; see Remark 2.6. 

In the sequel we will frequently use the bivariate random variable (A(t),A(s)). 
It obviously obeys a two-dimensional Normal distribution with zero mean 
and covariance matrix S(s,i). With T(s,t) := Cav(A(t), A(s)), this covari- 
ance matrix is given by 



Gaussian sources have the conceptual problem that the possibility of neg- 
ative traffic is not ruled out, as opposed to "classical" input processes, such 
as (compound) Poisson processes or on-off sources. However, in queueing 
theory a key role is played by functionals of the arrival process, which are 
well defined, regardless of whether the input stream corresponds to nonneg- 
ative traffic or not. Consider, for instance, the stationary distribution of a 
queue fed by a single source, emptied at rate c, given by the well-known 
formula sup t>0 (A(—t,0) — ct). Clearly, the distribution of such functionals 
can still be evaluated for Gaussian input; see, for example, Norros' pioneer- 
ing work for fBm [26], or [15]. We remark that such an approach leads to 





n — >oo 77, 
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nonnegative queue lengths in tandem systems with Gaussian inputs — this 
will follow directly from representation (10). For priority systems it is ex- 
plained in detail in [25], Section 2.3, how negative queue lengths can be 
avoided (a discrete-time version of the priority discipline is introduced, in 
which negative traffic can annihilate queued traffic). 

2.3. Sample-path large deviations. The analysis in the next sections re- 
lies on a sample-path large-deviations principle (LDP) for centered Gaussian 
processes. This section is devoted to a brief description of the main theorem 
in this field, (the generalized version of) Schilder's theorem [5]. However, we 
start by recalling (the multivariate version of) the well-known Cramer's the- 
orem; see [10], Theorem 2.2.30. We let (•, •) denote the usual inner product: 
(a, b) := a T b = Ya=i a %h- 

Theorem 2.1 (Multivariate Cramer). Let X{ G M d be i.i.d. d-dimensional 
random vectors, distributed as a random vector X with finite moment- generating 
function Ke^' X ' ( for all £ M. d ). Then the following LDP applies: 

(a) For any closed set F C M. d , 




(b) For any open set G C M d , 




where the large deviations rate function A(-) is given by 

(3) A(x) := sup((0,x) -logEe^' X) ). 

Remark 2.2. Consider the specific case that X has a multivariate Nor- 
mal distribution with mean vector \x and {d x d) nonsingular covariance 
matrix S. Using logEe< e ' X > = (0,/x) + ±6> T £6», it is not hard to derive that, 
with (x - fi) T = (xi — jLti, . .. ,x d - (jb d ), 

(4) r = S" 1 (x-^) and A(x) = \{x - fifZ' 1 ^ - p), 
where 9* optimizes (3); it is well known that A(-) is convex. 

We now sketch the framework of Schilder's sample-path LDP, as estab- 
lished in [5], see also [11]. We adopt the notation and setup of [1, 23]. 
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Consider the n i.i.d. centered Gaussian processes Ai(-), as introduced in 
Section 2.2. Define the path space £1 as 

n := {uj : R -> R, continuous, cj(0) = 0, lim f , , = lim ^ , . = ol, 

^ t-»oo 1 + |t| t^-oo 1 + |t| J 

which is a separable Banach space by imposing the norm 

II II 

||w||fi :=sup . 
teR 1 + Fl 

In [1] it is pointed out that -Aj(-) can be realized on Q under the assumption 
that 

(5) lim — - = for some a < 2. 

t^oo t a 

We assume assumption (5) to be in force throughout this paper. 

Next we introduce and define the reproducing kernel Hilbert space R C 
fi — see [3] for a more detailed account — with the property that its elements 
are roughly as smooth as the covariance function r(s, •). We start from a 
"smaller" space R*, defined by 

R* := |w :R-^ R, w(-) = 5^air(si,-),Oi, Si €R,neN|. 

The inner product on this space R* is, for Lu a ,uJb £ R* , defined as 

In n \ n n 

(6) (uj a ,uj b ) R := (J2air(si,-),j2 b j r ( s jr)) = a i b j T ( s ^ 

\i=l j=l I R i=lj=l 

notice that this implies {T(s, •), t))n = T(s, t). This inner product has 
the following useful property, which we refer to as the reproducing kernel 
property: 

n In \ 

(7) ^(t)=^a i r( Si ,t) = (^a l r( Si ,-),r(i,-)) =(u{-),T{t,-)) R . 

i=l \i=l I R 



From this we introduce the norm \\oj\\r := \/ (u,co)r. The closure of R* under 
this norm is defined as the space R. Now we can define the rate function of 
the sample-path LDP: 

^H^llfj, if oj E R, 
oo, otherwise. 

For the Gaussian sources introduced in Section 2.2, the following sample- 
path LDP holds. 



(8) 



Theorem 2.3 (Generalized Schilder). The following sample-path LDP 
applies: 
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(a) For any closed set F <z£l, 



1. _/l n 



limsup-logP - V^j(-) G F < - inf Ji 



(b) For any open set G C £1, 



Urn inf- log Pf - Y'M-) G G | > - inf 



A difficulty of Schilder's theorem is its "implicitness," as only in special 
cases the rate function /(•) can be explicitly minimized over the set of inter- 
est. The authors of [1] succeed in exploiting the reproducing kernel property 
to give a sample-path analysis of overflow in a single FIFO queue (with de- 
terministic service rate nc) fed by Gaussian inputs. With Q n denoting the 
stationary buffer content, they derive 

1 (b + ct) 2 

lim - log F(Q n > nb) = - inf V , / . 

rwoo n Wn ' t>o 2v(t) 

It is elementary to show the existence of a minimizing t. On the one hand, 
it holds that \mitio{b + ct) 2 /2v (t) = oo [due to v(0) = 0]. On the other hand, 
(5) implies that there is a (3 < 2 such that v(t) < tr eventually, and hence 
also limt_» ^ (6 + ct) 2 /2v(t) = oo. Notice that even in the FIFO setting, there 
is not necessarily uniqueness of the optimizing t; see, for instance, [21], [1], 
Section 3.7, or [20], Example 5.2. 

If t* denotes a minimizing t, the corresponding path is 

< 9 > ™ = -^po + ^ - - v{tt) ~ Zpj^ (» + ^ 

This path corresponds to a buffer that starts to fill at time — t* , and reaches 
overflow at time 0; it is not hard to check that /*(—£*) = — b — ct* , and 
/*(0) = 0, as desired. Notice that the path is in R (in fact even in R*). 
If there is a unique optimizing path in the target set (i.e., the set of all 
paths leading to overflow), it is usually referred to as the most likely path to 
overflow. It has the interpretation that, given that the rare event of overflow 
happens, with high probability it happens according to this trajectory. Also, 
t* has then the interpretation of the most likely duration of the busy period 
preceding overflow. (Notice that, in this FIFO setting, there is not necessarily 
uniqueness; see, e.g., Section 3.7 in [1] or Example 5.2 in [20], or, in a non- 
Gaussian setting, [21].) 
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2.4. Many-sources scaling. In this section we show that the probability 
of our interest can be written in terms of the "empirical mean process" 
n 1 J27=i The following lemma exploits the fact that we know both a 

representation of the first queue Qi <n (in steady state) and a representation 
of the total queue Q\, n + Q2,n (in steady state). Let to := bj (ci — c 2 ). 

Lemma 2.4. f(Q 2:n > nb) equals 

P|3t>t :VsG (0,t):- Vii(-f-s) >b + c 2 t-c 1 s\ 

\ n ti J 

Proof. Notice that a "reduction principle" applies: the total queue 
length is unchanged when the tandem network is replaced by its slowest 
queue; see [4, 14]. More formally: Q 1>n + Q 2 , n = sup t>0 (Y,i=i Ai(-t, 0) - 
nc 2 t). Consequently we can rewrite 

Q2.11 = (Ql.n + Q2,n) — Ql,n 

(10) 

= sup V^i(-t,0) - nc 2 t - sup Vi,(-s,0) - ncis . 
*>°\i=i / s >°\i=i J 

It was shown (see Lemma 5.1 of [31]) that the negative of the optimizing 
t (s) corresponds to the start of the last busy period of the total queue 
(the first queue) in which time is contained. Notice that a positive first 
queue induces a positive total queue, which immediately implies that we can 
restrict ourselves to s £ (0,t). Hence P(Q2,n > nb) equals 

P|3t>0:VsG(0,t):iV AA-t, -s) > b + c 2 t - as] . 

V n ti J 

Because for s f t the requirement 
1 n 

■ ^ t, -s) > b + c 2 t - c±s 



n . 



reads > b + (c 2 — c\)t, we can restrict ourselves to t > to- We can interpret 
to as the minimum time it takes to cause overflow in the second queue (notice 
that the maximum net input rate of the second queue in a tandem system 
is ci — c 2 ). □ 

The crucial implication of the above lemma is that for analyzing F(Q 2t n > 
nb), we only have to focus on the behavior of the empirical mean process. 
More concretely, 



(11) P(Q 2 , n >n6)=P 
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Fig. 1. Graphical representation of the overflow set. For different values oft, the curve 
b + c-it — c\ (t — s) has been drawn. Overflow occurs if there is at > to such that the empirical 
mean process lies, for s £ (0,t), above the corresponding curve. 

where the set of "overflow paths" S is given by 

S:={fen:3t>t ,\fse(0,t):f(-s)-f(-t)>b + c 2 t-c 1 s}. 



Remark 2.5. A straightforward time-shift shows that the probability 
that the empirical mean process is in S coincides with the probability that 
it is in T, with 

(12) T:={/GQ:3t>t ,VsG (0,t):/(s) > b + c 2 t - ci{t - s)}. 

However, the set T is somewhat easier to interpret, see Figure 1. For different 
values of t [i.e., t 2 > t\ > to = b/{c\ — c 2 )], the line b + c 2 t — c\{t — s) has been 
drawn. The empirical mean process n _1 Y^i=i ^i(') is i n T if there is a t > t$ 
such that for all s G (0, t) it stays above the line b + c 2 t — c\ (t — s) . Notice that 
T resembles the set corresponding to the probability of long busy periods 
in a single queue, as studied in [27]. 

Remark 2.6. As indicated above, our results are for centered sources, 
but they can be translated easily into results for noncentered sources. Then 
the traffic generated by Gaussian source i in the interval [s,t) is A(s,t) + 
/j,(t — s), where A(s,t) corresponds to a centered source; here < \i < min{ci, c 2 } 
and s < t. Let q(fi, ci,c 2 ) be the probability that the second queue exceeds 
nb, given that input rate \x and service rates c\ and c 2 are in force. From (10) 
it follows immediately that q(fi,ci,c 2 ) = q(0,c± — fi,c 2 — fi), and hence we 
can restrict ourselves to centered sources. 
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3. Analysis. In this section we analyze the logarithmic asymptotics of 
^(Q2,n > nb). In Section 3.1 we show that the decay rate in (2) exists, of 
which we derive a lower bound in Section 3.2. It turns out that this lower 
bound has an insightful interpretation, which is given in Section 3.3. Sec- 
tion 3.4 presents conditions under which the lower bound is tight (meaning 
that the decay rate and lower bound match). Finally, in Section 3.5 we prove 
and explain some properties of the most likely path that we found. 

3.1. Decay rate of the overflow probability. In this section we establish 
the existence of the decay rate (2) of P(Q2,n > Tib). We already saw in (11) 
that P(Q2 n > n b) can be rewritten as the probability that the empirical 
mean process is in S (which is an open subset of Q). The existence of the 
decay rate follows from Schilder's result (Theorem 2.3), by showing that S 
is an /-continuity set, that is, that the infima of /(•) over S and S match. 

Theorem 3.1. 

J=inf J(/) = inf /(/). 
fes fes 

The proof of Theorem 3.1 can be found in the Appendix. 

3.2. Lower bound on the decay rate. The main result of this section is a 
tractable lower bound on J, which is given in Theorem 3.2. Observe that 

S= (J p| S s ' f withcS 5 '* :={/€«: /(-a) -/(-*)> 6 + pzt-cis}. 

t>t se(o,t) 

Hence we are interested in the decay rate of the union of intersections. The 
decay rate of a union of events is simply the minimum of the decay rates of 
the individual events. The decay rate of an intersection is not standard. In 
the next theorem we find a straightforward lower bound on this decay rate. 
Define 

LI 8 '* : = {/ € n : -/(-*) > b + c 2 t; f(-s) - f(-t) >b + c 2 t- as}. 



Theorem 3.2. The following lower bound applies: 
(13) J> inf sup inf 1(f). 

*>*>.6(0,t) 



Proof. Clearly, 



J=inf inf /(/). 

t>to /-pO 5s, t 
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Now fix t and consider the inner infimum. If /(— s) — f(—t) > b + c 2 t — c\s 
for all s G (0, t), then also (/ is continuous) /(— s) — /(— t) > 6 + c 2 i — c%s for 
all s G [0, t] . Hence, 

P| cS s '*c p| w^c^ 1 
se(o,t) se[o,t] 

for all r G (0, t), and consequently 

inf /(/)> inf /(/). 

Now take the supremum over r in the right-hand side. □ 

Theorem 3.2 contains an infimum over / G Z^ s '*. In the next lemma we 
show how this infimum can be computed. Recalling (4), the bivariate large- 
deviations rate function of 

/A ^(-t,o) A M-t,-s) \ 

\ i=l n i=l n / 

is, for y,z eM. and t > 0, s G (0,t), given by A(y, z) := |(y, z)£(t — s,t) _1 (y, z) T . 
We also introduce the following quantity, which plays a key role in our anal- 
ysis: 

k(s,t) :=E(A(-s,0)\A(-t,0) = b + c 2 t) 

(14) 

= E(A(s)\A(t) = b + cat) = ^^(6 + c 2 t). 

Under the following assumption, the infimum over XA s,t can be simplified 
considerably. The same assumption will be useful when deriving tightness 
conditions in Section 3.4. 

Assumption 3.3. \A>(-) G C 2 ([0, oo)) is strictly increasing and strictly 
concave. 

Lemma 3.4. Under Assumption 3.3, for t > tQ and s G (0, t), 

;„ f T(f\_ T („ f \._f A ( b + c 2t,b + c 2 t-c 1 s), ifk(s,t)>c 1 s, 
f&t.,tV } L[S ' Z) - \{b + c 2 t) 2 /2v{t), tfk(s,t)< Cl s. 

PROOF. Observe that 



(15) 



> >o + c 2 t;> > o + c 2 t - cis . 



n t - ; n 

, i=i j=i 
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Hence we can use Theorem 2.1, yielding 

/ inf/(/)=infA(y,z), 

where the last infimum is over y >b + C2t and z > b + c^t — c\s. Using that 
A(-, •) is convex, this problem can be solved in a standard manner. It is easily 
verified that the contour of A that touches the line y = b + C2t does so at 
z-value 

T(t-s,t) n 
v{t) 

also the contour that touches z = b + C2t — c\s does so at y- value 

r(t-s,t),, 
yo ■= — 7z — r \ b + c 2* ~ c i s )- 

v(t — s) 

We first show that it cannot be that yo > b + C2t, as follows. If yo > b + C2t, 
then the optimum would be attained at (yo,b + C2t — c±s). Straightforward 
computations, however, show that y$ > b + C2t would imply that [use T(t,t — 
s)<Vv{t)v{t-s)} 

(16) (\/^)- ^Jv{t-s))(b + c 2 t) > ^t)c lS . 

This inequality is not fulfilled for s = (0 ^ 0) nor for s = t (b + C2t ^ cit for 
t > to)- As the left-hand side of (16) is convex (in s) due to Assumption 3.3, 
whereas the right-hand is linear (in s), there is no s € (0,i) for which the 
inequality holds. Conclude that yo > b + C2t can be ruled out. Two cases are 
left: 

(A) Suppose zq > b + C2t — c±s, or, equivalently, k(s,t) < c\s. Then (b + 
C2*, zq) is optimal (see the left panel of Figure 2), with rate function 
(b + C2t) 2 /2v(t), independent of s. 

(B) In the remaining case (where y® <b + C2t and zq < b + C2t — cis) the 
optimum is attained at (b + C2t, b + C2t — c\s), that is, the "corner point"; 
see the right panel in Figure 2. This happens if k(s,t) > cis, and gives 
the desired decay rate. 

This proves the stated. As an aside we mention that if k(s,t) = cis, then 
both regimes coincide: A(b + C2t, b + C2I — c\s) = (b + C2t) 2 /2v (t). □ 

Corollary 3.5. Under Assumption 3.3, the following lower bound ap- 
plies: 

J > inf sup T(s, t). 

i>io se(o,t) 
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z 




z 



b + c 2 t 

—as 



b + c 2 t 

-as 




b + at 



b + c 2 t 



y 



Fig. 2. Contour lines of the (two-dimensional) rate function; the objective function is to 
be minimized over the shaded region. 

3.3. Interpretation of the lower bound. The results of the previous sec- 
tion have a helpful interpretation, leading to two regimes for values of c±. 
For ci smaller than some critical link rate cf , we show in Corollary 3.7 that 
the lower bound of Corollary 3.5 can be simplified considerably. 

We start by drawing a parallel with the single-node FIFO result, as dis- 
played in (1). There, t has to be found such that 



is minimized. Let tf denote an optimizing argument t. L c (t) can be inter- 
preted as the cost of generating b + ct in an interval of length t, and t^ as 
the time duration yielding the "lowest cost." 

Now we turn to our tandem setting, and in particular to the result of 
Lemma 3.4. Computing the minimum of A(y, z) over its admissible region, 
we saw that, under Assumption 3.3, in both cases the optimizing y was equal 
to y = b + C2t. On the contrary, for the optimizing z there were two possible 
regimes. 

Now recall the representation (14) of k(s, t) as a conditional mean, and (15). 
The result in Lemma 3.4 essentially states that in the regime k(s, t) < c±s the 
most likely realization of J27=i 1,0) > nb+nc2t yields Ya=i Ai(—t, —s) > 
nb + nc2t — nc\s (with high probability, n large). In the other regime, k(s,t) > 
c\s, the most likely realization of J27=iAi(— 1,0) > nb + nc2t does not auto- 
matically yield Ya=i ^ ~ s ) <nb + nc2t — nc\s (with high probability, n 
large); fulfilling the second constraint in (15) requires additional "cost." 

The next decomposition result follows immediately from Lemma 3.4 and 
the above. 



Corollary 3.6. For s e (0,t), we have T(s,t) = L C2 (t) + L(s\t), with 

, . . _ m&x 2 {E(A{s)\A(t) = b + C2t)-c lS ,0} 
[S]): ~ 2Var(A(s)\A(t)=b + c 2 t) 



L c (t) : 



{b + ct-EA(t)) 2 
2VarA(t) 
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(17) 



max 2 {£;(s, t) — cis, 0} 



2Var (A(s)\A(t) =b + c 2 t)' 



Similarly to the interpretation of the single-node FIFO result, we can inter- 
pret T(s, t) as the cost of generating the required amount of traffic. Denoting 
by s* and t* optimizing arguments in Corollary 3.5, the intuition is as fol- 
lows: 

(A) "Cost component" L C2 (t) is needed to generate b + c 2 t in the interval 
(— t, 0]. By taking the infimum over t (to get £*) we find the most likely 
epoch to meet the constraint. 

(B) "Cost component" L(s\t) is required to make sure that no more than c\s 
is generated in the interval (— s,0], conditional on the event A{— i,0) = 
b + c 2 t. We can interpret s* as the epoch at which most effort has to be 
done to fulfill this requirement. This is of course reflected by the fact 
that in Corollary 3.5 we have to take the supremum over all s in (0,t). 
Evidently, if k(s,t) < C\s for all s E (0,i), this cost component is 0. 

For large values of ci, k(s,t) will be smaller than c±s for all s £ (0,t), since 
it does not depend on c\. As argued above, in this case the second term in 
Corollary 3.6 vanishes. If this holds for the t that maximizes the first term, 
that is, tc , then 



It implies that, for these large values of c%, the lower bound on J of Corol- 
lary 3.5 coincides with the result of a single-node FIFO queue with service 
rate c 2 . The intuition behind this is that essentially in this regime all traffic 
entering the first queue is served immediately, and goes directly into the 
second queue; traffic is not "reshaped" by the first queue. If c\ < cf , then 
the first queue does play a role in delaying and reshaping the traffic before 
entering the second queue, as we will see in the next section. 

Corollary 3.7. For all c\ > cf , (18) applies. 



(18) 



t>to se ( ,t) 



inf sup T(s,t) = L C2 {t^ 2 ). 



This clearly holds for all c\ larger than 



cf := inf{ci|VsG (0,t£) : fc(s,tf 2 ) <c x s\ 
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3.4. Tightness of the decay rate. Corollary 3.5 is a lower bound on the 
decay rate J. Of course, such a bound is only useful if it is relatively close 
to the actual decay rate, or, even better, coincides with it. In the latter case 
we say that the lower bound is tight 

In Section 3.2, we have derived a lower bound on J by replacing the decay 
rate of an intersection of events by the decay rate of the least likely of these. 
It is important to observe that if the optimum path in this least likely set 
happens to be in all the sets of the intersection, then the lower bound is 
tight. 

More specifically, let s* and t* be optimizers in the lower bound of Corol- 
lary 3.5. Clearly we can prove tightness of the lower bound by showing that 
the most probable path in IA S '* is in S (or S; use Theorem 3.1). In our 
analysis we distinguish between (A) c\ > cf , and (B) c\ < cf . 

Regime (A) [c\ larger than the critical service rate). In this situation, 
we know from Corollary 3.7 that the lower bound in Corollary 3.5 reduces 
to the decay rate in a single FIFO queue. The next result follows easily. 

Theorem 3.8. Under Assumption 3.3, if c\ > cf, then 

J= inf sup T(s,t) = L C2 (if ), 
*>*o s e(o,t) 

and a most probable path in S is 

(19) r(r) = -E(A(r,0)\M-C0) = b + c 2 t^ 2 ). 

Proof. As shown in Section 3.3, in this regime t* =t£, whereas the 
choice of s* is irrelevant [as c\ > cf implies L(s\t*) = for all s £ (0,t*)]. 
Notice that it is now sufficient to show that /* G S, or f*£S (use Theo- 
rem 3.1). We claim that /*(•) G S, or more precisely, that there exists t>tQ 
such that for all s G (0,t) it holds that f*(-s) - /*(-t) > b + c 2 t-c 1 s. This 
follows because, by definition of cf , for all s G (0,i*), 

/*(-«) " /*(-**) = nM~t\ s)\A(-t*, 0) = b + c 2 t*) 
= b + c 2 t* - k(s, t*)>b + c 2 t* - cis. 

This completes the proof. □ 

We want to stress that the above theorem holds for all Gaussian processes, 
regardless of the specific shape of the variance function. Consequently, the 
result is also valid for long-range dependent processes, such as fractional 
Brownian motion. 



16 M. MANDJES AND M. VAN UITERT 

Regime (B) (ci smaller than the critical service rate). We follow the 
same approach as in Regime (A): first we derive (in Lemma 3.10) a most 
probable path in 14 s '* , and then we verify (in Theorem 3.11) whether this 
path is in S. It turns out that we have to impose certain additional conditions 
to make the lower bound of Corollary 3.5 tight. We proceed by two technical 
lemmas; the proof of Lemma 3.9 is given in the Appendix. 

Lemma 3.9. Under Assumption 3.3, if c\ < cf , then k(s*,t*) > c\s*. 

Lemma 3.10. If k(s,t) >c±s, then a most probable path inlA s,t is 
(20) f(r) = -E(A(r, 0)\A{-t, 0)=b + c 2 t, A(-s, 0) = c lS ), 

with norm A(b + c 2 t, b + C2t — c\s). 

Proof. Using standard properties of conditional multivariate Normal 
random variables, we see that f(r) equals 

-et(s,t)T(-r,t)-8*(s,t)T(-r,s) 

We finish the proof by applying Lemma 3.4, and observing that 

Ill/Ill = T(a, t) = A(6 + oat, b + c 2 t- as), 
which is a matter of straightforward calculus. □ 

Before presenting our tightness result for the case c\ < cf, we introduce 
some new notation: 

(i) For r\,T2 < 0, 

EA(n,r 2 ) :=E(A(n,r 2 )\A(-t*,0) = b + c 2 t*), 

with Var(-) and Cov(-,-) defined similarly. Also, v(r\) := Var A(r\, 0) and 
f (n, r 2 ) := Cov(A(n , 0), A(r 2 , 0)). 

(ii) For r £ (— 1*,0) we define the functions 

EA(r,0) + Cl r m{r) F(r,-s*) 
m{r) := , m[r) := — -, p(r) :~ 



y/v(r) m(-s*)' ^v(r)v(-s*) 

Theorem 3.11. Suppose 

(22) m(-s)<p(-s) for alls G (0,t*). 

Under Assumption 3.3, if c\ < cf, then 

J=inf sup T(s,t) = A(6 + c 2 t*,6+ c 2 t* - as*), 
t>to se(0,t) 
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and a most probable path is 

/*(r) = —E(A(r, 0)\A(-t*,0) = b + c 2 t\ A(-s\ 0) = as*). 

Proof. As in Theorem 3.8, we have to show that /*(•) is in S. This is 
done as follows: 

/*(-*) " /(-**) = E(A(-t*, -s)\A{-t\ 0) = b + c 2 t\ A{-s*, 0) = cis*) 
= b + c 2 t* - E(A(-s, 0)\A(-s*, 0) = cis*) 

= b + c 2 t* - EA(-s,0) - r( ~ g '~f ) (cis* - EA(-a*,0)). 

!)(-S ) 

Now it is easily seen that (22) implies that /*(— s) — f(—t*) >b + c 2 t* — c\s 
for all s£ (0,t*). 

Due to Lemma 3.9, k{s* ,t*) > cis*. With Lemma 3.10, the expression for 
J follows. □ 

Although the condition (22), required in Theorem 3.11, is stated in terms 
of the model parameters, as well as known statistics of the arrival process, 
it could be a tedious task to verify it in a specific situation. The next lemma 
presents a somewhat more transparent necessary condition for (22). 

The intuition behind the lemma is the following. Observe that both p(-) and 
attain a maximum 1 at r = — s*. For p{ ) this follows from the observa- 
tion that p{r) is a correlation coefficient; for m(-) from Corollary 3.6 and 
Lemma 3.9. Then a necessary condition for (22) is that in s* the curve 
m(-) is "more concave" than p(-). The proof of the lemma is given in the 
Appendix. 

Lemma 3.12. A necessary condition for (22) is 

(23) m"(-s*)<p"(-s*), 
or equivalently, 

(24) 0£(s*,t*)(u"(t* -s*) - v"(s*)) + 6* 2 (s\t*)(v"(0) -v"(s*)) > 0. 

Condition (24) has an insightful interpretation, which will be given in the 
next section. 

3.5. Properties of the input rate path. So far, we have analyzed paths / 
of the cumulative amount of traffic injected into the system. In this section 
we turn our attention to the first derivative of /, which can be interpreted 
as the path of the input rate of the queueing system. As before, we have to 
consider two regimes: (A) c\ > cf , and (B) c\ < cf ; let Assumption 3.3 be 
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in force. Consider the paths /* as identified in Theorems 3.8 and 3.11, and, 
more specifically, their derivative <?*(•) := (/*)'(•)• I n case (A), with t* = tf 2 , 
and r G (-t*,0), 

9\r) = b -±^-tv'{r + t*) + v'{-r)\ 
whereas in case (B) it turns out that, with r G (— t*, — s*], 

and with r G [— s*,0), 

If v'(0) = 0, we show below that the path g*(-) has some nice properties. 
Notice that the requirement v'(0) = holds for many Gaussian processes. It 
is not valid for standard Brownian motion (Bm), since then v(t) = t, but the 
special structure of Bm allows an explicit analysis, see Section 4.1. Fractional 
Brownian motion (fBm), with v(t) = t 2H , has t/(0) = only for H G (^, 1]; 
see Section 4.2. 

Proposition 3.13. If ci > cf and v'(0) = 0, then g*(0) = g*(-t*) = c 2 . 
Proof. Notice that, due to (1), t* satisfies 

The stated follows immediately from v'(0) = 0. [As an aside, we mention 
that <?*(•) is symmetric in —t*/2.] □ 

Just as we exploited properties of t* in the proof of Proposition 3.13, we 
need conditions for s* and t* in the regime c\ < cf . These are derived in the 
next lemma. 

Lemma 3.14. If c\ < cf , then s* and t* satisfy the following equations: 
2c 2 = et(s*,e)v'(t*) + 0£(aV*)(«'(t*) - v'{f - «*)), 
2ci = 0%(8*,t*)i/(s*) + 9i(s*,t*)(v'(s*) + v'(t* - 8*)). 

Proof. By Lemma 3.9, k(s*,t*) > as*. Observe that T(s,t) = A(b + 
C2t, b + C2t — c\s) can be rewritten as 



(25) 9 T x(s, t) - i6» T S(s, t)8 where x(s, t) := 



b + c 2 t 

C\S 
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here we abbreviate 6 = (0*(s, t), 9 2 (s, t)) T . We write d t and d s for the partial 
derivatives with respect to t and s, respectively. The optimal s* and t* 
necessarily satisfy the first-order conditions, obtained by differentiating (25) 
to t and s, and equating them to 0. Direct calculations yield 

hc 2 \ _ ( m d t e 2 \ ... , ... / \elv'{t) + d t r( s , t)exh 



hci)~\d s e 1 d s e 2 j^ {s ' t)e x( - s ^ + \ieiv'( s ) + d s v(s,t)e 1 e 2 

The second equality in (21) provides x(s, t) = E(s, t)6. Now the stated follows 
directly. □ 

Proposition 3.15. // c\ < cf and v'(0) = 0, then (i) g*(—t*) = c 2 , 
and (ii) g*(—s*) = c\. Also, the necessary condition (24) is equivalent to 
(g*Y(-s*)>0. 



Proof. Claims (i) and (ii) follow directly from v'(0) = and Lemma 3.14. 
The last statement follows directly after some calculations. □ 



Proposition 3.15 can be interpreted as follows. The second queue starts 
a busy period at time —t*. During this trajectory, the first queue starts to 
fill at time — s* and is empty again at time 0, if the conditions of Theo- 
rem 3.11 apply. It is also easily seen that the necessary condition (24) has 
the appealing interpretation that (<?*)'(— s*) > 0: the input rate path should 
be increasing at time — s*. 



3.6. Some remarks. 



Remark 3.16. In our lower bound we replace the intersection over 
s G (0, t) by the least likely event of the intersection. Under condition (24) 
the occurrence of the least likely event implies all the other events in the 
intersection, with high probability [in the sense that /* E 14 s '* implies that 
/* e U s ' 1 * for all s G (0,i*)]. The examples in Section 4 show that (22) is 
met for many "standard" Gaussian models, but not always. If there is no 
tightness, a better lower bound can be obtained by approximating the in- 
tersection by more than just one event: 

J> inf sup inf /(/), 

*>*o se(o,t) m /ew s -* 

where s = (s±, . . . , s m ), and the "multiple-constraints set" U s,t is defined by 
:={/£fi: -f(-t) >b + c 2 t; 

f(-Si) - f(-t) >b + c 2 t- cisi, for i = 1,.. . ,m}. 



Obviously, the lower bound becomes tighter when increasing m. 
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Remark 3.17. The approach we have followed in this section to analyze 
the two-node tandem network can be easily applied to an m-node tandem 
network, with strictly decreasing service rates, that is, c\ > ■ ■ ■ > c m — nodes 
i for which a < Cj+i can be ignored, see [4, 14, 16]. Note that J2i=i Qi,n is 
equivalent to the FIFO queue in which the sources feed into a buffer that 
is emptied at rate c^. This means that we have the characteristics of both 
Yh=i Qi,n and YaLi Qi,n, which enables the analysis of Qm t n, just as in the 
two-node tandem case. 

4. Examples. One of the reasons for considering Gaussian input pro- 
cesses is that they cover a broad range of correlation structures. Choosing 
the variance function appropriately, we can make the input process exhibit- 
ing, for instance, long-range dependent behavior. In this section we do the 
computations for various variance functions. We also discuss in detail the 
condition in Theorem 3.11. 

4.1. Standard Brownian motion. The variance function for Brownian 
motion (Bm) is given by v (t) = t. Using (1), it is easily found that tf 2 = b/c 2 . 
According to Corollary 3.7, cf is the largest value of c\ such that for all 

*e(o,t£)> 

^(b + c 2 t^)- Cl s<0, 

that is, cf = 2c2- Hence, using Theorem 3.8, we have for c\ > 2c 2 that J = 
2bc 2: with a constant input rate g*(r) = 2c 2 for r G (— ,0) and g*(r) = 
elsewhere. 

Now we turn to the case where c\ < 2c 2 - The optimizing s* and t* are 
determined by solving the first-order equations for s and t; see Theorem 3.11. 
We immediately obtain that t* = b/ {c\ — c 2 ) and s* = 0. Obviously, for this 
regime the service rate of the first queue does play a role. The most probable 
input rate path reads g*(r) = c\, for r £ (— i*, 0) and g*(r) = elsewhere. It 
is easily verified that the most probable path /*(•) is in S, making the decay 
rate as found in Theorem 3.11 tight. In other words, 

be 2 

J = A(b + c 2 t\ b + c 2 t* - cis*) = — i — r. 

2(ci - c 2 ) 

Observe that, interestingly, Bm apparently changes its rate instantaneously, 
as reflected by the most likely input rate path. This is a consequence of the 
independence of the increments. 
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4.2. Fractional Brownian motion. The variance function for fractional 
Brownian motion (fBm) is given by v(t) = t 2H , where H is the so-called 
Hurst parameter. For H > i this corresponds to long-range dependent traf- 
fic. Now (1) gives 



t F 



b H 



By Theorem 3.8, 



J=- 



1/ 6 \ 2 " M / C2 -^ M 



2\l-H V # 



for all ci > Cj . Unfortunately, for general H there does not exist a closed- 
form expression for cf. Now turn to the case c\ < cf . Lemma 3.12 states 
that (24) is a necessary condition for tightness to hold. Observe that v"(t) = 
(2H- l)2Ht 2H - 2 and hence v"(0) = oo. It is easily checked that 0£(s*,t*) < 
0, which implies that in this case (24) is not satisfied. Therefore the lower 
bound on J is not tight. 



4.3. M/G/oo input. A versatile traffic model is the so-called M/G/oo 
input process. In this model sessions arrive according to a Poisson process 
with rate A, and stay in the system for some random duration D. During this 
period they generate traffic at a unit rate. By choosing specific session-length 
distributions D, both short-range and long-range dependent inputs can be 
modeled. For more results on queues with M/G/oo input traffic processes, 
see, for example, [12, 30]. Below we approximate the M/G/oo inputs by 
their "Gaussian counterpart," that is, Gaussian sources with the same mean 
and variance as the M/G/oo input; this procedure is extensively motivated 
in [1, 2]. 

Let the mean session-length be finite, say 5, such that the mean input 
rate equals X5. We denote by the distribution function of D and 

by -Fi)r(-) the distribution function of the residual session-length, that is, 
FD r (x) = <5 _1 Jq(1 — Ff)(y))dy. We denote the corresponding densities by 
fn(-) and f D r(-). 

Let B(t) denote the amount of traffic generated by a single M/G/oo input 
in an interval of length t. We now show how to compute the variance v(-) 
of B(t). We will do this by first deriving the moment-generating function 
of B(t). In fact two types of sources contribute: 

1. Sources that were already present at the start of the interval. The number 
of these sources has a Poisson distribution with mean XS. Their residual 
duration has density fjjr ■(•); with probability (1 — Fjjr(t)) they transmit 
traffic during the entire interval. 
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2. Sources that arrive during the interval. Their number has a Poisson(At) 
distribution. Given that the number of these arrivals is k G NU {0}, their 
arrival epochs are i.i.d. random variables, uniformly over the interval 
(with density i" 1 ). Their duration has density /d(-). 

Straightforward computations now yield (cf. [22]) 

logE(e ei? W) = X5{M t {9) - 1) + \t(N t (6) - 1) 

with 

M t {9) := f e ex f D r{x) dx + e e \l - F D r(t)) 
Jo 

and 

N t (6):= C ( t -e^ x - u) f D {x-u)dxdu+ ( -e e ^ u \l - F D {t - u)) du. 

JO Ju t Jo t 

Taking the second derivative of the log moment-generating function (with 
respect to 9) and then substituting for 9, gives the variance v(t) of B(t): 

X5 (J X 2 f D r (X) dx+t 2 (l~ F D r (t))] 

+ x([ [ (x - u) 2 f D (x - u) dx du + [ (t-u) 2 (l- F D (t-u))du). 

\Jo Ju Jo J 

For fBm we could a priori rule out tightness of the lower bound due to 
v"(Q) = oo; see Lemma 3.12. For M/G/oo inputs we show in the following 
lemma that v"(0) is finite, even for heavy-tailed D. It implies that condi- 
tion (22) needs to be checked to verify tightness. 

Lemma 4.1. For 5 < oo and finite /d(-), both (i) v'(0) = and (ii) v"(0) < oo. 

Proof. Using standard rules for differentiation of integrals, 

v' (t) = XS2t(l - F D r(t)) + X f 2{t-u)(l-F D {t-u))du 

Jo 

and hence v'(0) = 0. Similarly, 

V"(t) = 2X6(1 -F D r(t)-tf D r(t)) 

+ 2X f (1 — Fd (t — u) — (t — u)f D (t - it)) du 
Jo 

POO 

= 2XJ (l-F D (s))ds. 

Hence, v"(0) = 2X5 < oo. □ 

Now we consider some examples of session-length distributions. In all the 
examples we take b = 0.5, A = 0.125, 5 = 2 and C2 = 1. 
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Fig. 3. M/exp/oo input process. 



Exponential. Using the above formula for v(-), we get 



v(t) = 2X5 3 - - 1 + exp 
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I — g*(0 

1.2 - 
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-5 -4.5 -4 -3.5 -3 -2 5 -2 -1.5 -1 -0 5 

*- time 



Fig. 4. Input rate path for M/exp /oo input process. 

Notice that v (•) tends to a straight line for large t (corresponding to short- 
range dependence). Numerical computations then give cf = 1.195. Tak- 
ing ci = 1.1 results in s* = 4.756, t* = 5.169 and m(r), p(r) as given in 
Figure 3. The upper panel of Figure 3 shows m(r) and p{r) for r G (— 1*,0), 
whereas the lower panel magnifies the graph around — s* . We see that in- 
deed m(-) < p(-) on the desired interval, so the decay rate is tight. A corre- 
sponding input rate path is given in Figure 4, which satisfies the properties 
as indicated in Proposition 3.15. 

Hyperexponential. In case D has a hyper exponential distribution, with 
probability pi G (0, 1) it behaves as an exponential random variable with 
mean , with i = 1,2 and p\ + p2 = 1- It is easily verified that 

v(t) = 2\^(v 1 t - 1 + e-^) + 2X^ 2 t - 1 + e"^), 

with V2 =P2/{$ —px/v\). As in the exponential case, v(-) is asymptotically 
linear. For p\ = 0.25 and v\ = 5, we find cf = 1.173, and s* = 4.700, t* = 
5.210, when using c\ = 1.1. Also for this example m(-) < p(-), and hence there 
is tightness; the graph looks similar to Figure 3. A corresponding input rate 
path is given in Figure 5. 
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-5 -45 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 

— » time 

Fig. 5. Input rate path for M/H2/oo input process. 



i I l r I l r I l r 




-6 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 

«. time 



Fig. 6. Input rate path for M/Par/oo input process. 
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Pareto. If D has a Pareto distribution, then F(D >t) = (l/(l + i)) Q . The 
variance function is given by 

v ^ = 7^ u 9 2A vi - (* + 1 ) 3_Q + ( 3 " 

(3 — a)(2 — a)(l — a) 

with a= (1 + 8)/6, excluding 5 = 1 or i. Notice that we have a = 1^, yield- 
ing v(i) ~ t\/t, which corresponds to long-range dependent traffic. Numeri- 
cal calculations show that cf = 1.115, and for c\ = 1.1 we obtain s* = 4.373, 
t* = 5.432. Again m(-) is majorized by An input rate path is given 
in Figure 6. We empirically found that there is not always tightness in the 
M/Par/oo case. If b is larger, for instance 6 = 1, then (22) is not met. 



5. Priority queues. In Section 3 we analyzed overflow in the second 
queue of a tandem system. This analysis was enabled by the fact that we 
had explicit knowledge of both the first queue and the total queue. In the 
present section we use the same type of arguments to solve the (two-queue) 
priority system. 



5.1. Analysis. We consider a priority system with a link of capacity nc, 
fed by traffic of two classes, each with its own queue. Traffic of class 1 does 
not "see" class 2 at all, and consequently we know how the high-priority 
queue Qh,n behaves. Also, due to the work-conserving property of the sys- 
tem, the total queue length Qh,n + Qe,n can be characterized. Now we are 
able, applying the same arguments as for the tandem queue, to analyze the 
decay rate of the probability of exceeding some buffer threshold in the low- 
priority queue. This similarity between tandem and priority systems has 
been observed before; see, for instance, [13]. 

We let the system be fed by n i.i.d. high-priority (hp) sources, and an 
equal number of i.i.d. low-priority (lp) sources; both classes are indepen- 
dent. We assume that both hp and lp sources are Gaussian, and satisfy the 
requirements imposed in Section 2. Define the means by \ih and //^, and the 
variance functions by Vh(-) and vi(-), respectively; also \i := fj,h + l^e (where 
fx < c) and v(-) := ityi(-) +vt{-). We note that in this priority setting we can- 
not restrict ourselves to centered processes. We denote the amount of traffic 
from the ith hp source in (s,t], with s < t, by Ah : i(s,t); we define A^ n (s,t) 
analogously. Also Th(s,t),T£(s,t) and Rh,Re are defined as before. 



Remark 5.1. Notice that this setting also covers the case that the num- 
ber of sources of both classes are not equal. Assume, for instance, that there 
are na lp sources. Multiplying m and vt{-) by a and applying the fact that 
the Normal distribution is infinitely divisible, we arrive at n i.i.d. sources. 
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In the tandem situation we could, without loss of generality, center the 
Gaussian sources. It can be checked easily that such a reduction prop- 
erty does not hold in the priority setting, since there is no counterpart 
of Remark 2.6. Hence we cannot assume without loss of generality that 

Analogously to Lemma 2.4, we obtain that P(Q£,n > nb) equals 

/ 1 n 1 n \ 

P 3t>0:Vs>0:-V A hji (-t, -s) + - V A^-t, 0) > b + c(t - s) . 

Let J p be the exponential decay rate of P(Q^ n > nb); analogously to Theo- 
rem 3.1 it can be shown that this decay rate exists. Similarly to the tandem 
case, with /(•) = (//>(•), f e (-)), 

S'/ :={/eflxfi: f h (-s) - f h {-t) - M-t) > 6 + c(t - s)}, 
V ; p V fh{s)- fh(-t)- ft(-t)>b + c(t-s) J 



/ 4 i=l *=1 / t>0s>0 / 



Theorem 5.2. T/ie following lower bound applies: 
(27) J p > inf sup inf /(/), 

<>° s>0 fell*'* 

with f h (t) := f h {t) - fi h t, fe(t) := f e (t) - u. e t and 

i(f) ■■= III Alii + k\\ftf Rl - 

The infimum over / £ Wi 3 '* can be computed explicitly, as in Lemma 3.4. 
As the analysis is analogous to the tandem case, but the expressions are 
more complicated, we only sketch the procedure. Again there is a regime in 
which one of the two constraints is redundant. Define 

k p {s,t) := E(A h (s)\A h {t) + A e (t) = b + ct). 

Using the convexity of the large-deviations rate function, it can be shown 
that, if 

E(A h (t -s) + A e {t)\A h {t) + A e (t) = b + ct)>b + c{t - s), 

only the first constraint in (26) is tightly met; it is equivalent to require that 
k p (s,t) < cs. [If k p (s,t) > as, either both constraints in (26) are met with 
equality, or only the second constraint is met with equality; exact conditions 
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for these two cases are easy to derive, but these are not relevant in this 
discussion.] As before, under k p (s,t) < cs, we obtain the decay rate 

(28) M<fW - g±fezffg , 

compare the FIFO queue with link rate nc; in the other cases the expressions 
are somewhat more involved. Denote by t F the value of t > that minimizes 
the right-hand side of (28). 

Similarly to the tandem case, there is a regime (i.e., a set of values of the 
link rate c) in which J p coincides with the decay rate of an FIFO queue. 
In this regime, which we call regime (A), conditional on a large value of 
the total queue length, it is likely that the hp queue is empty, such that all 
traffic that is still in the system is in the lp queue. Hence, for all c in 

(29) {c\Vs>0:k p (s,t F ) <cs} 
we conclude 

p ~ 2v{t F ) 

If c is not in the set (29), we can use the methodology of Section 3 to find a 
condition under which the lower bound of Theorem 5.2 is tight; we call this 
regime (B). 

Remark 5.3. In the tandem case, we found that the FIFO result holds 
for ci > cf , whereas it does not hold for c\ < cf; the threshold value cf was 
found explicitly in Section 3.3. In the priority setting there is not such a 
clear dichotomy. Consider, for instance, the situation in which both types 
of sources correspond to Brownian motions; Vh{t) = X^t, vi{t) = X(t and 
A := Xh + A^. Define 



'/^ + ^(c-/Uft) 2 - 



Then straightforward calculus yields that for (A/j — A^)c < XhifJ-h + 2/^) — 
XtHh, regime (A) applies (i.e., the FIFO result holds): 

_ 26(c-/i) 
Jp- A , 

whereas otherwise we are in regime (B): 

_ 6(5 - fit) _ 
P ~ X, ' 

this is shown by verifying that the lower bound of Theorem 5.2 is tight for 
the specific case of Brownian motion input. Using fi^ + fi£ < c, it can be 
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verified easily that this implies that for A^ < Xe the FIFO solution applies, 
whereas for A/i > A^ only for 

c< Afc(/i fe + 2fi e ) - X e n h 
Aft - A^ 

the FIFO solution applies. 

5.2. Discussion. Large deviations for priority queues have been studied 
in several papers. We mention here the work by Mannersalo and Norros [23] 
and Wischik [33]. We briefly review their results, and compare them with 
our analysis. Our lower bound then reads 

jjP := inf supT„(s, t) with T p (s, t) : = inf 1(f)- 

Just as we did, Mannersalo and Norros [23] identify two cases. They get the 
same solution for our regime (A), that is, the situation in which, given a 
long total queue length, the hp queue is relatively short; see also Berger and 
Whitt's [6] empty buffer approximation. 

In regime (B) the hp queue tends to be large, given that the total queue 
is long. To prevent this from happening, [23] proposes a heuristic that min- 
imizes /(/) over 

(30) {/£Oxfi:3t>0: -fh(-t) - fe(-t) >b + ct; -f h (-t) < ct}. 

Because regime (B) applies, the optimum paths in the set (30) are such that 
the constraints on / are tightly met; consequently (30) is a subset of Up*. 

Hence the resulting decay rate, which we denote by Jp l \ yields a lower 
bound, but our lower bound will be closer to the real decay rate: 




Remark 5.4. In the simulation experiments performed in [23], the lower 

bound Jp(b) is usually close to the exact value. Our numerical experiments 
(cf. the examples on the tandem queue in Section 4) show that the hp buffer 
usually starts to fill shortly after the total queue starts its busy period. This 
means that in many cases the error made by taking s = t is relatively small. 
It explains why the heuristic based on set (30) performs well. 

Wischik [33] focuses on discrete time, and allows more general traffic than 
just Gaussian sources. Translated into continuous time, in regime (B), his 

lower bound on the decay rate Jp U ^ (Theorem 14) minimizes 1(f) over 
{/G^xO:3t>0:3s>0: -fa(-t) - ft(-t) >b + ct; 

(31) 

-fh(s)<cs}; 
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again a straightforward comparison gives that our lower bound Jp is closer 
to the actual decay rate: 

4 m) :=inf inf TJs,t) < inf sup T p (s,t) = J® . 

Remark 5.5. Recent work by Mannersalo and Norros [24] suggests that 
a similar approach could work for a queue operating under the Generalized 
Processor Sharing (GPS) scheduling discipline. For two classes of traffic 
(both with n sources), sharing a resource with link capacity nc and two 
buffers, the model is parametrized by the weights </>i,</>2 6 [0,1], summing 
to 1. If both queues are nonempty, both classes receive their guaranteed ser- 
vice rates ncfiic and nfac, respectively. If one class does not use all its band- 
width, it can be taken over by the other class in a work-conserving manner. 
For more details on the system mechanics for GPS, see, for instance, [28, 29]. 

Consider the probability that the first queue exceeds level nb, under the 
assumption that the mean input rates of both classes are smaller than their 
respective guaranteed service rates. Notice that the backlog of type 2 does 
not exceed that of an FIFO queue with link rate nfoc- This suggests that, 
in self-evident notation, the decay rate is well approximated by the infimum 
of 1(f) over 

{/ G x n :3t > 0:Vs > :-/«(-*) - f {2) {-t) - f {2) (-s) 

>b + ct — C(J)2S}. 

Reasoning heuristically (see also [24]), it is not likely that (i) queue 2 is 
nonempty at the start of the busy period preceding overflow of queue 1, 
(ii) there is traffic left in queue 2 at the epoch queue 1 reaches overflow. 
This would lead to a minimization over 

{/eQx ft:3t>0:Vs>0:-/ (1) (-i) - f {2) {-t) >b + ct; 

-f {2) (s)<cfos}, 

compare the sets U s,t (as identified for the tandem system) and lAt^ 1 (priority 
system). A lower bound for this decay rate is again found by taking the 
infimum over t > and the supremum over s > 0, as before. 

APPENDIX 

Proof of Theorem 3.1. To prove Theorem 3.1, we first present an 
auxiliary result (cf. [27], Proposition 4.2). The set S is open; we now deter- 
mine its closure. Define 

S l := {/ G : Vs G (0,t) : f(-s) - f(-t) >b + c 2 t- cys}, 

S s ^ :={fen: f(-s) - f(-t) >b + c 2 t- cis}. 
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Lemma A.l. The closures of S l and S are characterized as follows: 

(32) 5* = {/ G : Vs G (0, t) : f(-s) - f(-t) >b + c 2 t- as}, 

(33) S = [j 5*. 

t>t 

PROOF. We first prove (32). "C" is obvious: 

se(o,t) »e(o,t) 

Now consider "D." Let / be in the right-hand side of (32). Define, with 
y + :=max{0,y} and y~ :=min{0, y}, 

f n (u):=f(u) + -(u- +t)+. 
n 

It is easy to see that (i) ||/ — f n \\n — ► 0, and (ii) /„ G 5*; here (ii) follows 
from 

/n(-«) " /«(-*) = /(-«) " /(-*) + ^(t - s) > 6 + c 2 t - c lS 

for s G (0, i). This proves (32). 

Next we show (33). Again we establish two inclusions. "5" is done by 
picking an arbitrary / from the right-hand side: 

(i) Suppose there is a t > to such that / G S l ; then we can reuse the 
above argument: take an / from the right-hand side of (33), and show that 
there is a sequence f n in S such that ||/ — f n \\n 0. This is exactly as 
before. 

(ii) Suppose / is only in the union in the right-hand side of (33) for 
t = to; then we have to show that / can be approximated by an f n G S to+Sn , 
with ||/ — / n ||n — ► 0, and 5 n := 1/n. This is done by the following sequence: 

|7(t), fori>-t , 
fn(t) = \ cx(t + t ) + f(-to), for te[-to-S n ,-t ], 

{ -ci5 n + f{t) + /(-t ) - /(-to - S n ), for t<-t - S n . 

Now f n £S to+s ", as can be seen as follows. For s G (0,to), using / G 5'° in 
conjunction with (32), 

fn(s) ~ fn(~to ~ S n ) = f(s) + Ci5 n - /(-to) 

> b + C2to - cis + ci5 n >b + C2(to + 5 n ) - c\s, 
due to ci > C2- For s G [to, to + S n ), similarly, 

fn(-s) - fn(~to ~ S n ) = -C X S + Cit + C\5 n 

= -cis + b + c 2 t + c\8 n >b + c 2 (t + S n ) - c\s. 



•32 



M. MANDJES AND M. VAN UITERT 



Now concentrate on "C"; take /g5: 

(i) Hence there is a sequence 5 such that ||/ — f n \\n 0. Because 
f n G 5, there is a sequence of epochs t n (all of them strictly larger than to) 
such that f n G S tn . 

(ii) t n is bounded. This can be seen as follows. Clearly, due to f n G S tn , 

fn{-s) - fn{-t n ) >b + C 2 t n - CiS 

for all s G (0, t n ). Hence also — f n (—tn) > b+C2t n (let s J, 0), and consequently 

11/ /II ^ fi-tn) ~ fn(-t n ) ^b + C 2 t n f{~t n ) 

117 " fJn - TTT n - T+tn~ + T+1~' 

Suppose t n were not bounded; letting n — > oo would lead to a contradiction: 
> C2 [use that f(u)/(l + u) — > for u— >• oo]. 

(hi) Hence we can pick a subsequence t nk such that t nk goes to some 
finite limit > to for k — > oo. Now / G 5*°°, since, for all s G (0, t^) and 
sufficiently large, 

/»*(-«) - fn k {-tn k ) >b + C 2 t nk - C X S. 

This proves the lemma. □ 

We now prove the theorem. Clearly, from Schilder's result, 

inf /(/)<•/ <inf :/(/). 
/es /e5 

To show the stated, we prove that the infima over S and 5 coincide. 

Let T be defined by (12), and let T* and T s '* be defined analogously to 
S l and 5 s '*; their closures are determined as in Lemma A.l. It is evident 
that the infima over S and T coincide. As mentioned above, our aim is to 
prove that the infima over S and S match, but it turns out to be more 
convenient to show that the infima over T and T match. 

This is done by choosing / from T n R arbitrarily, and showing that we 
can approximate it by a path in T . Clearly / G T** for some t* > to- Let 
C be an arbitrary path in R that is strictly positive in (0, i*], and define 
fn '■= f + C/ n - Then there is a t n > t* such that f n G T tn . This can be seen 
as follows. 

(i) First observe that f(s)>b + C2t* — c\(t* — s) for all s in the closed 
interval [0,t*]. Hence for all s G (0,t*], and n G N, 

fn{s) = f(s) + -C(S) > b + C 2 t* - Cl (t* - S). 

n 

(ii) As this inequality also holds for s = i*, we conclude that there is a 
t n > with f n (s) > b + C2t n — ci(t n — s) for all s G (0,t n ), or, equivalently, 
that f n G T*" . 
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Now notice that, for n — > oo, 



fn\\ 2 R 



f + -C 

n 



R 



which proves Theorem 3.1. □ 



2 

Ri 



Proof of Lemma 3.9. The lemma is proven in three steps. Notice 
that, as we are in regime (B), it holds that c\ < cf = K(tf ), with 



(34) 



rJ ~, \ k(s,t) 
K(t):= sup — ^-A 

se(o,t) s 



(i) In [9], Lemma 3.1, it is shown that, under (5) and Assumption 3.3, 
L C2 (t) is decreasing for t <t^ , and increasing for t>t^ . 

(ii) We now prove by contradiction that K(t*) > c\. Suppose K(t*) < c\. 
Then, by (34), for all s G (0,t*) it holds that k(s,t*) < c\s, and hence also 
sup sg( - 0jt *) L(s\t*) = 0; see (17). Now consider the decomposition of Corol- 
lary 3.6: 



(35) 



sup T(s,t) = L C2 (t) + sup L(s\t); 
sG(0,t) se(o,t) 



t* is minimizer of this expression. Because k(s,t) is continuous, also for 
the closed interval [0,t*] it holds that k(s,t*) < c±s. Hence it is possible to 
decrease t* such that the first term in the right-hand side of (35) decreases 
(as we approach t^ 2 from above, see Step 1), while the second remains 0. 
Hence the sum of both terms decreases, implying that t* cannot be optimal. 
So it cannot be that both 



K(t*) < ci and t* > t 



< ci implies t* = t 



Similarly K(t*) < c\ rules out t* <t^ 2 . Hence K(t* 
However, K(t^) > c\. Contradiction. 

(hi) Notice that K(t*) > c±, in conjunction with (17), directly implies 
that k(s*,t*) > c\s*. This proves the lemma. □ 



Proof of Lemma 3.12. First we show that (23) holds. As noted earlier, 
both m(-) and p(-) have a maximum 1 at — s*. This means that (23) is 
necessary to enforce m(r) < p(r) for r in a neighborhood of —s*. 

Next we show that (23) is equivalent to (24). First multiply both m(-) 
and p(-) by h(-), where 

h(r) :=J^^(EA(-s\0)-c lS n. 
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Since h(r) : (— t* , 0) — > R+, this yields the requirement ir(r) < n(r) for all 
r G (-i*,0), with 

ffr — s*l 

vr(r) :=EA(r,0) + cir and n(r) := -±r r-(EA(-s*, 0) - as*). 

v(—s*) 

Recall that m(-) and p(-) have the same function value and derivative at 
— s*. It is easy to derive that this implies that (m ■ h)(—s*) = (p ■ h)(—s*) 
and (m • h)'(—s*) = (p ■ h)'(—s*). Therefore, the necessary condition be- 
comes 7r"(— s*) < n"(—s*). 

Using standard formulas for conditional means of multivariate Normal 
random variables, 

v(t*) 



leading to 



2 



^L(EA(r,0) + Cl r) 



b + c 2 t* 



Assuming r < — s*, 

= . _ ^ _ v{-r) + v{s*) - v{-r - s*) _ T(-r,t*)T(s* ,t*) 
(r ' S j " 2 «(**) 

such that 



u"(s*) - v"(0) v"(s*) - v"(t* - s*) v(s*, t*) 



2 2 v(t*) 

It can be checked that the same result holds when the derivative is calculated 
for r > — s*. Now it is a straightforward but tedious computation to prove 
that this implies that ir"(—s*) < n"(—s*) is equivalent to (24). □ 

Acknowledgment. The authors thank Krzysztof Dehicki (CWI) for valu- 
able discussions. 



REFERENCES 

[1] Addie, R., Mannersalo, P. and NORROS, I. (2002). Most probable paths and per- 
formance formulae for buffers with Gaussian input traffic. European Transactions 
on Telecommunications 13 183-196. 

[2] Addie, R., Neame, T. and Zukerman, M. (2002). Performance evaluation of a 
queue fed by a Poisson Pareto burst process. Computer Networks 40 377-397. 

[3] Adler, R. (1990). An Introduction to Continuity, Extrema, and Related Topics for 
General Gaussian Processes. IMS, Hayward, CA. MR1088478 

[4] Avi-Itzhak, B. (1965). A sequence of service stations with arbitrary input and reg- 
ular service times. Management Science 11 565-571. MR174382 



TANDEM AND PRIORITY QUEUES 



35 



Bahadur, R. and Zabell, S. (1979). Large deviations of the sample mean in general 

vector spaces. Ann. Probab. 7 587-621. MR537209 
Berger, A. and Whitt, W. (1998). Effective bandwidths with priorities. 

IEEE/ACM Transactions on Networking 6 447-460. 
Botvich, D. and Duffield, N. (1995). Large deviations, the shape of the loss curve, 

and economies of scale in large multiplexers. Queueing Systems 20 293-320. 

MR1356879 

Courcoubetis, C. and Weber, R. (1996). Buffer overflow asymptotics for a buffer 
handling many traffic sources. J. Appl. Probab. 33 886-903. MR1401482 

Debicki, K. and Mandjes, M. (2003). Exact overflow asymptotics for queues with 
many Gaussian inputs. J. Appl. Probab. 40 704-720. MR1993262 

Dembo, A. and Zeitouni, O. (1998). Large Deviations Techniques and Applications, 
2nd ed. Springer, New York. MR1619036 

Deuschel, J.-D. and Stroock, D. (1989). Large Deviations. Academic Press, Lon- 
don. MR997938 

Duffield, N. (1998). Queueing at large resources driven by long-tailed M/G/oo- 
modulated processes. Queueing Systems 28 245-266. MR1628426 

Elwalid, A. and Mitra, D. (1995). Analysis, approximations and admission control 
of a multi-service multiplexing system with priorities. In Proceedings of the IEEE 
Infocom 463-472. 

Friedman, H. (1965). Reduction methods for tandem queueing systems. Oper. Res. 
13 121-131. MR175196 

HijSLER, J. and Piterbarg, V. (1999). Extremes of a certain class of Gaussian 
processes. Stochastic Process. Appl. 83 257-271. MR1708208 

Kella, O. and Whitt, W. (1992). A tandem fluid network with Levy input. In 
Queues and Related Models (I. Basawa and U. Bhat, eds.) 112-128. Oxford 
Univ. Press, London. MR1210563 

Kilkki, K. (1999). Differentiated Services for the Internet. Macmillan Technical Pub- 
lishing, Indianapolis. 

Kilpi, J. and NORROS, I. (2002). Testing the Gaussian approximation of aggregate 
traffic. In Proceedings of the Internet Measurement Workshop. 

Leland, W., Taqqu, M., Willinger, W. and Wilson, D. (1994). On the self- 
similar nature of Ethernet traffic. IEEE/ACM Transactions on Networking 2 
1-15. 

Mandjes, M. (2004). A note on the benefits of buffering. Stochastic Models 20 43-54. 
MR2036295 

Mandjes, M. and Kim, J. H. (2001). An analysis of the phase transition phenomenon 
in packet networks. Adv. in Appl. Probab. 33 260-280. MR1825326 

Mandjes, M., Saniee, I. and Stolyar, A. (2000). Load characterization, overload 
prediction, and load anomaly detection for voice over IP traffic. In Proceedings 
of the 38th Allerton Conference 567-576. 

Mannersalo, P. and NORROS, I. (2001). Approximate formulae for Gaussian pri- 
ority queues. In Proceedings of the ITC 17 991-1002. 

Mannersalo, P. and Norros, I. (2002). GPS schedulers and Gaussian traffic. In 
Proceedings of the IEEE Infocom 1660-1667. 

Mannersalo, P. and Norros, I. (2002). A most probable path approach to queue- 
ing systems with general Gaussian input. Computer Networks 40 399-412. 

NORROS, I. (1994). A storage model with self-similar input. Queueing Systems 16 
387-396. MR1278465 



36 



M. MANDJES AND M. VAN UITERT 



[27] NORROS, I. (1999). Busy periods of fractional Brownian storage: A large deviations 

approach. Advances in Performance Analysis 2 1-20. 
[28] Parekh, A. and Gallager, R. (1993). A generalized processor sharing approach to 

flow control in integrated services networks: The single node case. IEEE/ACM 

Transactions on Networking 1 344-357. 
[29] Parekh, A. and Gallager, R. (1994). A generalized processor sharing approach to 

flow control in integrated services networks: The multiple node case. IEEE/ACM 

Transactions on Networking 2 137-150. 
[30] Parulekar, M. and Makowski, A. (1997). M/G/oo input processes: A versatile 

class of models for network traffic. In Proceedings of the IEEE Infocom 419-426. 
[31] Reich, E. (1958). On the integro-differential equation of Takacs. I. Ann. Math. 

Statist. 29 563-570. MR93834 
[32] Weiss, A. (1986). A new technique for analyzing large traffic systems. Adv. in Appl. 

Probab. 18 506-532. MR840105 
[33] WlSCHlK, D. (2001). Sample path large deviations for queues with many inputs. Ann. 

Appl. Probab. 11 379-404. MR1843050 



CWI 

P.O. Box 94079 
1090 GB Amsterdam 
and 

Faculty of Electrical 
Engineering, Mathematics 
and Computer Science 

University of Twente 

P.O. Box 217 

7500 AE Enshede 

The Netherlands 

e-mail: michcl@cwi.nl 



CWI 

P.O. Box 94079 
1090 GB Amsterdam 
and 

Faculty of Economics 

and Business Administration 
Vrije Universiteit 
De Boelelaan 1105 
1081 HV Amsterdam 
The Netherlands 
E-MAIL: miranda@cwi.nl 



